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Abstract. The isotropic Dunkl oscillator model in the plane is investigated. 
The model is defined by a Hamiltonian constructed from the combination of 
two independent parabosonic oscillators. The system is superintegrable and its 
symmetry generators are obtained by the Schwinger construction using parabosonic 
creation/annihilation operators. The algebra generated by the constants of motion, 
which we term the Schwinger-Dunkl algebra, is an extension of the Lie algebra u(2) 
with involutions. The system admits separation of variables in both Cartesian and 
polar coordinates. The separated wavefunctions are respectively expressed in terms of 
generalized Hermite polynomials and products of Jacobi and Laguerre polynomials. 
Moreover, the so-called Jacobi-Dunkl polynomials appear as eigenfunctions of the 
symmetry operator responsible for the separation of variables in polar coordinates. 
The expansion coefficients between the Cartesian and polar bases (overlap coefficients) 
are given as linear combinations of dual —1 Halm polynomials. The connection with 
the Clebsch-Gordan problem of the sZ_i(2) algebra is explained. 
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1. Introduction 

This series of two papers is concerned with the analysis of the isotropic Dunkl 
oscillator model in the plane. The system will be shown to be superintegrable and 
the representations of its symmetry algebra will be related to different families of — 1 
orthogonal polynomials EJ S21 El El El 05] • 

A quantum system defined by a Hamiltonian H in d dimensions is maximally 
superintegrable if it admits 2d — 1 algebraically independent symmetry operators S{, 
1 ^ i ^ 2d — 1, that commute with the Hamiltonian 



where one of the operators is the Hamiltonian itself, e.g. S\ = "H. For a superintegrable 
system described by a Hamiltonian of the form 



where A is the Laplace-Beltrami operator, the symmetries Si will be differential 
operators. In this case, the system is said to be superintegrable of order I if I is 
the maximum order of the symmetry generators S{ (other than %). One of the most 
important features of superintegrable models is that they can be exactly solved. 

When i — 1, the constants of motion form a Lie algebra. When i = 2, the symmetry 
algebra is quadratic [TTj [T2~j [T3| [MJ E]. Substantial work has been done on these systems 
which are now well understood and classified (see [H S3 HQ 171 HQ US I2UI I2H E01 and 
references therein). Further developments in the study of integrable systems include 
progress in the classification of superintegrable systems with higher order symmetry 
|22l 1381 [39] . the examination of discrete/finite superintegrable models [25] and the 
exploration of systems involving reflection operators [TJ HH [151 El 123 [291 EH E21 EU [35] . 

We here examine the Dunkl oscillator in the plane. This model is possibly the 
simplest 2D system described by a Hamiltonian involving reflections and corresponds 
to the combination of two independent parabosonic oscillators. As will be shown, 
this system possesses many interesting properties. It is second-order superintegrable. 
Its symmetry algebra, which we term the Schwinger-Dunkl algebra, is obtained using 
parabosonic creation/annihilation operators in a way that parallels the Schwinger su(2) 
realization in the case of the ordinary 2-dimensional isotropic oscillator; the Schwinger- 
Dunkl algebra is an extension of the Lie algebra u(2) with involutions. The system 
admits separation of variables in both Cartesian and polar coordinates and its separated 
wavefunctions can be obtained explicitly in terms of the generalized Hermite, Jacobi and 
Laguerre polynomials. Furthermore, the study of this model and of the representations 
of its symmetry algebra will show remarkable occurrences of —1 orthogonal polynomials 
(OPs) families. The present paper is concerned with the exact solutions of the model, its 
superintegrability and the calculation of the overlap coefficients between the Cartesian 
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and polar bases. The second paper of the series will focus on the representations of the 
symmetry algebra and the connections with —1 OPs. 

Here is the outline of the paper. In Section 2, we define the Hamiltonian 
of the Dunkl oscillator and obtain its separated wavefunctions in Cartesian and 
polar coordinates. We also show that the symmetry operator responsible for the 
separation of variables in polar coordinates has the so-called Jacobi-Dunkl polynomials 
as eigenfunctions. In Section 3, we obtain the symmetry algebra of the model in terms of 
the parabosonic creation/annihilation operators. In Section 4, we show that the overlap 
coefficients between the Cartesian and polar bases are given by linear combinations of 
the dual —1 Hahn polynomials. In section 5, we exhibit the relationship between the 
Dunkl oscillator model and the Clebsch-Gordan problem of s/_i(2) [9| I40|. 

2. The model and exact solutions 

The isotropic Dunkl oscillator model in the plane is defined by the Hamiltonian 

n = -\ [{V^f + [V^f] + \[x 2 + y 2 ] , (1) 
where the operator is the Dunkl derivative 

Vfr=d Xi + ^Q-R Xi ), Xie{x,y}, (2) 
with I the identity operator and d x% = -Jr-. The operator R Xi is the reflection operator 

x i 

with respect to the plane Xi = 0. Hence the reflections in ([T]) have the action 

Rxf(x, y) = f(-x, y), R y f(x, y) = f(x, -y). 

In connection with the nomenclature of the standard harmonic oscillator, the model 
is called isotropic because the quadratic potential is 5*0(2) invariant. For the full 
Hamiltonian (CQ) to have this symmetry requires of course /j, x = fi y . Expanding the 
square of the Dunkl derivative, one finds 

{V ^ r = d 2 2 ^g t^ {I _ R ^ 

The Schrodinger equation 

= £V, (3) 

is manifestly separable in Cartesian coordinates. As shall be seen, even in the presence 
of reflections, (j3J) also admits separation in polar coordinates. Separation of variables 
in more than one coordinate system is a signal of superintegrability. This occurs for 
the Dunkl oscillator because reflections can be viewed as rotations. We provide below 
the exact separated solutions of (j3J). Note that when fi x = fi y = 0, the Hamiltonian (flj) 
corresponds to the standard quantum Harmonic oscillator in the plane. 
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2.1. Solutions in Cartesian coordinates 
Since the Hamiltonian (0Q) has the form 

W = *Hx + Hy, 

where % x is the Hamiltonian of the one-dimensional Dunkl oscillator, it is obvious that 
the solutions to ([3]) in Cartesian coordinates will be given by 

V(x,y)=Tl>(x)il>(y), £ = £ x + £y, 

where ip(xi) is an eigenfunction of the ID Hamiltonian with energy eigenvalue £ Xi . For 
the ID oscillator T-L x , the Schrodinger equation reads 

tfj"(x) + ^ il>\x) + (2S X - x 2 )tfj(x) - - R x )^{x) = 0. (4) 

X X A 

Since [H x , R x ] = 0, the eigenf unctions ip(x) may be chosen to have a definite parity 
R x ip(x) = s x ip(x) with s x = ±1. 

When s x = +1, we have R x ip + (x) = ip + (x) and the equation (TjJ has for (admissible) 
solutions 

I'M = J r( +?\ U9 , e-^L^\x*), 
]] T(n + n x + 1/2) 

where L^*\x) are the Laguerre polynomials [23J and where T(x) denotes the gamma 
function. The eigenvalues are given by 

£, = 271 + ^ + 1/2, ne {0,1,2,...}. 

When s x = —1, we have R x ip~(x) = —ip~(x) and the solutions to (T4|) are then 



y r(m + fx x + 3/2) 

with eigenvalues 

£ x = 2m + l + fix + 1/2, m e {0,1,2,...}. 



From the orthogonality relation of the Laguerre polynomials (lA.ip . it is easily seen that 
for fi x > —1/2, the eigenfunction i/;^ obey 

V>» 0*0 Vl>mY \x\ 2>lx dx = 5 nm 8 SxS , x , 

where x* denotes complex conjugation. From the above considerations, it is clear that 
the eigenstates of TL X can be labeled by a single integer n x whose parity is that of the 
corresponding wavefunction. For this purpose, we introduce the generalized Hermite 
polynomials [TJ |32j [33] 



"W'l = '-')V r (B+p+ ^i / ,) ''^W 
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where p = 0, 1. With this definition, the eigenf unctions of H x can be expressed as 

iP nx (x) = e- x2 / 2 HX : (x), n x G N, 

with energy eigenvalues £ x — n x + \i x + 1/2. 

The eigenf unctions of the one-dimensional Dunkl oscillator are thus normalized and 
orthogonal on the weighted L 2 space endowed with the scalar product 

poo 

(g\f)= / g*(x)f(x)\x\ 2 ^dx. (5) 



It is directly checked (see Appendix B) that the Dunkl derivative (j2J) is ant i- Hermit ian 
with respect to the scalar product (j3J). This establishes that the Hamiltonian (JTJ) is 
Hermitian. 

Using the above results for the one-dimensional Dunkl oscillator, it follows that the 
eigenstates of the full Hamiltonian (CQ) in the Cartesian basis satisfy 

W\n x ,n y ) = £\n x ,n y ), £ = n x + n y + \x x + p y + 1, (6) 

where n x , n y are non- negative integers. The wavefunctions have the expression 

* nx , ny (x,y) = e-^+y 2 y 2 H% : (x)H%(y), 

and they satisfy the orthonormality condition 

*n x ,n y (x,y)K> x ,n> y {x , y)\x\ 2 ^\y\ 2 ^ dx dy = 5 nx< 5 nv< , 

provided that \i x > —1/2 and n y > —1/2. For the ID case see also [26, 32J. 

2.2. Solutions in polar coordinates 
In the polar coordinate system 

x = p cos <f, V = P sin <p, 

the Hamiltonian (CTJ can be written as 



OO POO 



oo J — OO 



where A p has the expression 

A=- 1 - 
and where B v is given by 



d 2 p + - p d p 



- ~(Px+ Py)d p +l-p 2 , 
P ^ 
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For separation of the Dunkl Laplacian in higher dimensions see [31 El E]- The actions of 
the reflection operators are easily seen to be 

Rxf(p, <p) = f(p, 7T - <p), R y f(p, tp) = f(p, -tp). 

Upon substitution of the separated wavefunction *&(p, (p) = P(p)<&(<p) in ([3]), one obtains 
the pair of equations 

2 

Try 

A p P(p) - SP(p) + — 2 P{p) = 0, (8a) 



2p 2 



2 



B v *(<p) - — <%) = 0, (8b) 

where m 2 /2 is the separation constant. 

We start by examining the angular equation (j8b[) : it has the explicit form 



$" - 2U X tan p-fx y cot <pW - ^L_^l$ _ + m 2$ = . (9) 

cos"' </? sin 

Since ["H, = ["H, = 0, we shall label the eigenstates by the eigenvalues s x , s y = ±1 
of the reflection operators R x and R y . 

When s x = s y = +1, the equation Q has the (admissible) solution 

<£++ _ / (2n+fj, x +fi y )r(n+/i x +ri y )n\ p( fe -l/2, Hy-1/2) ( \ 
y 2r(n+^+l/2)r(n+n, + l/2) r n W> 

with x = —cos and where Pn l, ^\x) denotes the Jacobi polynomials [23] . This 
solution corresponds to the eigenvalue m 2 = 4n(n + \i x + /^) with n G N. 
When s x = s y = — 1, the solutions reads 



— y 2T(n+ Ata; +l/2)r(n+ AtH +l/2) S111 V 9 COS V 3 r n-l W> 

with variable x = — cos2<£> and eigenvalue m 2 = An(n + p, x + p y ), n6N. It is understood 
that P { ^\x) = and hence that $ _ = 0. 

When s x = +1 and = —1, the solution to equation Q) is given by 



<f,+- _ / (2n+f, x +r ly )r(n+r lx +^ y +l/2)(n~l/2)\ . p O*-l/2, ^+1/2) , x 

*n - V 2r( n +^)r(n+ ft +i) ^V-'n-i/a W> 

with variable x = — cos 2<£>, eigenvalue m 2 = 4n(n + p x + p y ) and where n takes only 
positive half-integer values n G {1/2, 3/2, 5/2, . . .}. 

Lastly, when s x = — 1 and s y = 1, the solution to the angular equation has the 
expression 



& -+ _ / (2n+^+ My )r(n+ M:c + My + l/2)(n-l/2)! _ p (/^+l/2, My -l/2) , x 
~ \l 2r(B+^+l)r{n +( ,„) C0S( ^n-l/2 W> 

with variable x = — cos 2</?, eigenvalue m 2 = 4n(n + /i x + /i y ) and where n takes only 
positive half-integer values n G {1/2, 3/2, 5/2, . . .}. 
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From the orthogonality relation of the Jacobi polynomials (IA.2j) . it is directly seen 
that the wavefunctions obey the orthogonality relation 

/ $^(<^)$™ s %) | cos<^| 2 ^| sin^| 2 ^ = 8 nrn 8 SxS ,J SySy . 
J o 

As seen from the above considerations, the value of the separation constant is always 
m 2 = 4n(n + p x + p y )- When the product s x s y = +1 is positive, n is a non-negative 
integer. When the product s x s y = — 1 is negative, n is a positive half- integer. 
We now examine the radial equation (IHaj) . It reads 

P'Xp) + ~ (1 + 2p x + 2p y )P\p) + Us - p 2 - P(p) = 0. 

This equation has for solutions 

niPj - y r(fc+2n+^+^+l) e ^ ^ >> 

with the energy eigenvalues 

£ = 2(k + n) + p x + fiy + 1, fceN. 

Using the orthogonality relation of the Laguerre polynomials, one finds that the radial 
wavefunction obeys 

POO 

/ P k (p)P k ,(p)p l+2 ^ydp = 5 kk , 
Jo 

Hence the eigenstates of the Hamiltonian ([T|) in the polar basis can be denoted 
\k, n; s x , s y ) and satisfy 

H\k,n;s x ,s y ) = £ \k,n;s x ,s y ), £ = 2(k + n) + p x + p y + 1, (10) 

where k e N is a non-negative integer and where n is a non-negative integer whenever 
the product s x s y = +1 is positive and a positive half-integer whenever the product 
s x s y = — 1 is negative. 

From the equations and (ITUj) . it is seen the states with a given energy 
£ = N + p x + p y + l exhibit aiV + 1-fold degeneracy. Here are the first few eigenstates : 



£ 


\Tl X , Tly) 


| k, Tt, S X j Sy) 


£ Q = 1 + P X + Py 
£ X = 2 + P X + Py 

£ 2 = 3 + P X + Py 
£ 3 = 4 + P X + Py 
£4 = 5 + P X + Py 


|0,0) 
|1,0),]0,1) 
|2,0>, H,l>, |0,2> 
|3,0), 1 2,1) , 1 1,2) , |0,3) 
4,0) [3,1), |2,2), 1 1,3) , 0,4) 


|0,0;++) 
|0,l/2;+-), |0,l/2;-+) 

1,0;++), |0,1;++), |0,1; ) 

l,l/2;+-),|l,l/2;-+),|0,3/2;+-),|0,3/2;-+) 
|2,0;++), |1,1;++), 1 1,1; ), |0,2;++), |0,2; ) 



The presence of these degeneracies can be attributed to the existence of a symmetry 
algebra that will be identified in Section 3. 
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2.3. Separation of variables and Jacobi-Dunkl polynomials 

As is seen from ( I8b[) . the separation of variables of the Schrodinger equation in polar 
coordinates is equivalent to the diagonalization of the operator B^. We thus have the 
following eigenvalue equation: 

2 

m 

B v ,\k,n;s x ,s y ) = —\k,n;s x ,s y ), m 2 = 4n(n + fi x + fi y ) , (11) 

where n G N when s x s y = 1 and n € {1/2, 3/2, . . .} when s x s y = —1. We shall consider 
the operator 

J 2 = i(xV%y-yV%*), 

which in polar coordinates reads 

J 2 = * [d v + H y cot (p(l- Ry) - /i x tan (p(I- R x )} . 

A simple computation shows that the square of the operator J 2 is related to B v in the 
following way: 

Ji = 2B^ + 2fi x/ jy(I-R x Ry). (12) 

Instead of the eigenvalue equation ( 1TTT) . we shall consider the one corresponding to the 
diagonalization of J 2 . 

J 2 F e (ip) = X e F e (<p) } (13) 

where e = s x s y = ±1; this extra label on the eigenvalues X e is allowed since R x Ry 
commutes with J 2 . It follows from (TTT]) and ( ITS]) that the square of the eigenvalues X e 
are given by 

\ 2 + = 4n(n + fi x + fiy), X 2 _ = 4(n + fi x )(n + fiy), (14) 

where n e N when e = 1 and n = {1/2, 3/2, . . .} when e = —1. Moreover, since s x s y = e, 
we have R x = eR y . To solve (I13j) . we consider the decomposition 

F e (<p) = m<p) + fr(<p), (is) 

where R y f^((p) = ±f^((p). It is directly seen that given the decomposition (fT5l) . the 
eigenvalue equation f fT3|) is equivalent to the system of differential equations 

d v [/+ + fr] + 2% cot <p f~ - ii x tan <p[(l- e)/+ + (1 + ej/f] = -iK [/+ + /"], 
0„ [-/+ + fr] + 2» y cot <p f: - fi x tan <p [(e - 1)/+ + (1 + e)/f] = -*A e [/+ - /f], 

where the second equation was obtained from the first one by applying These 
equations are easily seen to be equivalent to 

<Wr + 2 Py cot ¥>/f - tan <p(l + e)/ e " = -iXJ+, 
d vf? ~ tan (1 - e)/+ = -iX e f~. 
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2.3.1. The case e = 1 
When e = +1, one has 

d v f+ + 2fi y cot ipf+ - 2fi x tan <pf~ = -i\+f+ (16a) 
d*fZ = ~iKf+, (16b) 

Substituting (116bD in ( 1 16 aft yields the equation 

<9j/+ + (2 My cot up - 2fi x tan + X 2 + f+ = 0. 

Since = 4n(n + /i x + /^), we directly obtain the result 

with x = — cos2<y9, eigenvalues A + = ±2^Jn(n + fx x + /x y ) and neN. Consequently, for 
e = +, the eigensolutions of (IT2"]) are given by 

= ^(*) + ^d v P^(x), (17) 
where the eigenvalues are given by 



X + = ±2Jn(n + /i x + fi y ), fiGN, 

and where the parameters are a — fi x — 1/2, /3 — fi y — 1/2 and x = — cos2<^. When 
e = +1, R x = R y and the operator — %J 2 can be written as 

where 

A^ x ^ y = 2 2 ^ + ^)(sin|^|) 2 ^(cos^) 2 ^, 

with ^'(v 9 ) = ^^(v 9 )- This directly establishes that the polynomials defined by (TTT)) 
correspond to the so-called Jacobi-Dunkl polynomials studied in [2]. 

It is possible to express the eigenfunctions of ^ in terms of the wavefunctions, 
which are eigenfunctions of B v . By taking the derivative of equation Q with respect to 
(p for s x = s y = 1 and adjusting the normalization, one obtains 



9 V K + (y?) = 2Jn(n + n x + fi y )$ n (uj). 



Upon substituting this result in (fTTj) . one finds that for e = +1, the eigenfunctions 
F+((p) of J2 and their corresponding eigenvalues are given by 



F + (uj) = $++(<^) ± 1 <!>--&), A + = ±2 Jn(n + /i x . + fi y ). (18) 
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J8.5.& Thee = -1 case 

When e = — 1, the equations (12. 3p and (I2.3P become 

d v fZ + 1n y cot (f fZ = -iX-fZ, 
d v fZ - 2/j x tan <p /+ = -iX-fZ. 

The first equation is easily rewritten as 

<9j/+ + (2^ cot p - 2 Mx tan ^)^/+ + (A 2 _ - 4 Ma ^)/+ - = 0. 

Given the value of A^ defined in (TT4j) . we directly find 

fZ = cos V P^ /2 ^- 1/2 \x), fZ = i- (9,/+ - 2^tan^+) 

with x = — cos29?. For e = — 1, the eigenfunctions of J7"2 and their corresponding 
eigenvalues thus take the form 



F^cp) = f+(<p) ± fZ(cp), \- = ±2 y J(n + fj lx )(n + fj, y ), 

where x = — cos 2ip and where n is a positive half integer. In terms of the wavefunctions, 
a straightforward computation leads to the expression 



F_(y?) = $;+(v?) =F i A_ = ±2y> + ^(n + (19) 

Thus we have obtained the eigenfunctions of the operator J% in terms of the 
wavefunctions, which are the eigenfunctions of B v . 

3. Superintegrability 

In this Section we show that the Dunkl oscillator model in the plane is superintegrable. 
We recover the spectrum of the Hamiltonian algebraically using the parabosonic 
creation/annihilation operators and obtain the symmetries using the Schwinger 
construction. 

3.1. Dynamical algebra and spectrum 

We first consider the dynamical algebra of the Dunkl oscillator model. We introduce 
two commuting sets of parabosonic creation/annihilation operators 

K = ^ (xi + ) , Al = -L ( Xi - ) , 

where Xi G {x,y}. These operators have the non-zero commutation relations 
[A x , A x ] = I + 2fi x R x , [A y , A y ] = I + 2fi y R y . 
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In terms of creation/annihilation operators, the Hamiltonians 7-L x and T-i y have the 
expression 

U x = ^{A x ,Al}, n y = ^{A y ,Al}, (20) 

where {x,y} = xy + yx denotes the anti-commutator. Thus the 2-dimensional 
Hamiltonian of the Dunkl oscillator ([1]) has the simple form 

n= 1 -{A x ,Al}+ 1 -{A y ,Al}. 

In the preceding Section, the eigenvalues £ of % have been obtained analytically by 
solving the Schrodinger equation. They can also be obtained algebraically. Indeed, we 
have the additional commutation relations 

[U Xi ,A x% ]=-A Xi , [7^,4,1 = 4, (21a) 
{A Xi ,R Xi } = {Al,R Xi } = 0, [U Xi ,R Xi ] = 0, (21b) 

where Xi G {x,y}. It is easily seen from the relations ( 1201) . ( 12 lap and (121bl) that the 
operators H Xi , A Xi , A x . and R Xi realize two independent copies of the parabosonic 
algebra which we have related to sZ_i(2) in [40j. It follows directly from the above 
commutation relations that 

£x = n x + fx x + 1/2, S y = n y + fi y + 1/2, n x ,n y eN. 

A direct computation shows that the action of the ladder operators A x , A\ on the 
Cartesian eigenbasis \n x ,n y ) is given by 

Al\n x ,n y ) = \J[n x + l]^\n x + l,n y ), A x \n x ,n y ) = \J[n x ]^ x \n x - l,n y ), (22) 
and that of the reflection R x by 

Rx\^X1^"y) ( 1) \ r Tlxi'l r l"y) i 

where denotes the 'mu- numbers': 

[n] M = n + A*(l -(-!)")• 

Analogous formulas hold for the action of A y , A y and R y . 

As noted previously, the spectrum of the Hamiltonian TL has the form 

E N = N + /jl x + Hy + 1, N EN, 

and exhibits & N + 1-fold 'accidental' degeneracy at level N. These degeneracies will 
be explained in terms of the irreducible representations of the symmetry algebra of the 
Dunkl oscillator. 
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3.2. Superintegrability and the Schwinger- Dunkl algebra 

We now exhibit the symmetries of the Hamiltonian (CQ). Let us consider the operator 

J* = \{a x , 4} - \{Ay, 4} = \ (n x - Hy) . 

It is clear that [H, J3] = and that J3 is the symmetry corresponding to separation 
of variables in Cartesian coordinates. Following the Schwinger construction |36j, we 
further introduce 

J2 = Yi (AlA y - A X A\) . 

A direct computation shows that J 2 is also a symmetry, i.e. [H, J 2 ] = 0. In addition, 
expressing the operator J 2 in terms of Dunkl derivatives shows that 

J 2 = l-( x V^-yV^), 

and hence J 2 = —J 2 j2\ it is thus seen from (1121) that J 2 is associated to the separation 
of variables in polar coordinates. To obtain the complete symmetry algebra, we define 
a third operator which also commutes with %: 

Ji = \ (AlA y + A X A\) . 

A direct computation show that the symmetry operators of the Dunkl oscillator in the 
plane satisfy the following algebra 

{J u R Xi } = 0, {J 2 ,R Xl } = 0, [J 3 ,R Xl ] = 0, 

[J2,J3]=iJi, [Ja,Ji]=iJ2, 

[Jl, J2) =i[J3 + Ji [y-xRx + HyRy) ~ U{yL x R x ~ fX y R y )/2] , 

with R x — R 2 y = I; %i ^ {x, y} and where the Hamiltonian % is a central element. We 
shall refer to the algebra generated by J±, J 2 , J3, R x , R y and % as the Schwinger-Dunkl 
algebra sd(2); special cases of it have appeared in other contexts jH [H]. It is easily seen 
that sd(2) is a deformation of the Lie algebra u(2) by the two involutions R x , R y . The 
Schwinger-Dunkl algebra admits the Casimir operator [9] 

C = J-y + J 2 + t/3 4" —fX X R X + —fXyRy + fX X fXyR X Ry, 

which commutes with all the generators. A direct computation shows that in the present 
realization, the Casimir operator C takes the value 

4 4 

Since H is a central element, we can define 



C = C - H 2 /4 + 1/4, 
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and thus C = in this realization. 

The irreducible representations of the Schwinger-Dunkl algebra sd(2) can be used 
to account for the degeneracies of the Hamiltonian (CQ). We shall postpone this study 
for the second paper of the present series. Note that upon taking \x x = \x y = in the 
Schwinger-Dunkl algebra, the involutions cease to play an essential role and one recovers 
the well-known su(2) symmetry algebra of the standard quantum harmonic oscillator in 
the plane. 

4. Overlap Coefficients 

In this section, we obtain the expansion (overlap) coefficients between the Cartesian 
and polar bases. These expansion coefficients are denoted by (k,n; s x , s y \n x ,n y ) . It is 
clear that the coefficients will vanish unless the involved states \k, n; s x , s y ) and \n x , n y ) 
belong to the same energy eigenspace. The states in the polar basis are the eigenstates 
of the operator given in and satisfy 



with n a non-negative integer whenever the product s x s y = 1 and a positive half-integer 
otherwise. We can consider the relation 



and expand the action of B v on the Cartesian basis to obtain a recursion relation for 
the overlap coefficients. It will prove more convenient to investigate first the overlap 
coefficients between the Cartesian basis and the eigenbasis of a new operator Q related 
to J2. The eigenstates of this new operator Q will then be expanded in terms of the 
polar basis \k, n; s x , s y ) to obtain the desired result. For this part, it is convenient to 
separate the two eigenvalue sectors corresponding to the value of the product s x s y = ±1. 

4-1. Overlap coefficients for s x s y = +1 

We start by expressing the energy eigenstates in the polar basis in terms of the 
eigenstates of J2. As is seen from (118j) . the eigenvectors of the operator J2 with 
eigenvalues that we denote \n, ++) j 2 and \n, H — ) j 2 are given by 



Btp\k] n] s x , Sy/ — Tnl^j ^ 



;s x ,Sy), 7„ = 2n(n + fj, x + 



/ 7i(^5 s y \n x ^ n y ^ (k^ Ti, s x ^ s. 



v 




K 



ri 



+ 



n 





for n 7^ 0. For n 



0, one has 



0,++) j2 = |M;++) 



.+ _ 
— 



0. 
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We also recall that R y \n, ++) j 2 = \n, H — ) j 2 . We now introduce the operator Q defined 
by 

Q = iJlRx — H>xRy — l^yRx — (^/2)R x R y . (23) 

The relevance of the operator Q will become clear in Section 5 when the connection 
between the Schwinger-Dunkl algebra and the Clebsch-Gordan problem of sl-±(2) will 
be established. In the sector s x s y = +1, we have R x = R y and Q may be written as 

Q = lJ 2 Ry - fi x Ry - VyRy - (1/2)1 

For n^O, the eigenvalues and eigenvectors \n, +±)q of Q are found to be 
\n, ++)s = ^ {(n\n, ++)j 2 + \n, +-)j 2 ) , g+ = -2n - /i x - n v - 1/2, 



and 



\n, +-)q = -j= (-( n \n, ++)j 2 + \n, , q n = 2n + fi x + fi y - 1/2, 



where we have defined 



Cn 



Hx + Hy - 2iy/n(n + fi x + Jiy) 
2n + fi x + jj,y 



This amounts to the diagonalization of a 2 x 2 matrix. We note that ( n (* = 1. When 
n — 0, one has directly 

|0, ++)s = |0, ++) j 2 , qt = -fi x -fi y - 1/2. 

It is possible to regroup the eigenvalues of Q into a single expression. We have 

q e = (-l) i+1 (£ + fi x + lu y + l/2), 

and the eigenvectors are given by 

Mq = ^={Q \h + (i - S j0 )\j, +-)j 2 ), 

\Q2j+l)Q = -J= ( - 0+1 \j + h++)j2 + \ j + 1 ^+-)j 2 )- 

In the previous formulas, it should be understood that for the vector \q ) the 
normalization factor \/2 is not needed. 

Having introduced the operator Q, we examine the overlap coefficients between its 
eigenstates and the eigenstates of H in the Cartesian basis for a given energy level Sn- 
In the sector s x s y — +1, the possible levels take the energy values 



£ jv = N + n x + \i y + 1, 
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where iV is an even integer. The eigenspace Sn is spanned by the vectors 

|0, N), |1, N - 1), . . . , |m, N - m), . . . , \N, 0). 
We shall denote the overlap coefficients by 

(q e \m,N-m) = M^ N . 
To obtain the expression for the expansion coefficients N , we start from the relation 

q,M^ N = (q e \Q\m,N-m). (24) 
In terms of the parabosonic creation/annihilation operators, the operator Q reads 

Q = {A X A\ - A\A y )R x - (fx x + ^ y )R x - (1/2)1 (25) 
Upon substituting (|2"5]) in f)24p and using the actions (122 p . there comes 
Qi m L,n = An+i M m+1>JV + B m M f mtN + A m M f m _ 1N , 

where 

A m = (-lr^im^iN-m + l]^, B rn = (-l) rn+1 (fi x + fi y ) - 1/2. 

It follows that the overlap coefficients N can be expressed in terms of polynomials 
'P m {qi)- Indeed, if we define 

ML, N = Ml N v m (qt), 

with Vo(qg) = 1, it transpires that V m {qi) are polynomials of degree m in the variable 
qi obeying the three-term recurrence relation 

qiV m {qi) = A m+1 Tm+i{qi) + B m V m {qi) + A m V m -i(qt). (26) 
Upon introducing the monic polynomials V m (qi): 

V m (qe) = - A t-, 

the recurrence relation (1261) becomes 



Qe V m (q e ) = V m+ i(qi) + B m V m (qi) + U m V m -x{qe), (27) 

where 

U n = A 2 n = [m}, x [N-m + l}, y . (28) 

Comparing the formulas ( T27|) and (|28|) with the formula (IA.3[) of Appendix A, it is 
seen that the polynomials V m (qi) correspond to the monic dual —1 Hahn polynomials 

Q n (xi; a, (3; N). We thus have 

V m {q,) = 2- m Q m {x, 1 a 1 fi-N) 
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where the parameter identification is given by 

a = 2/iy + N + 1, P = 2fi x + N+1, 

and the variable X£ takes the values 

x e = (-l) m (2£ +2^ + 2^ + 1), £ = 0,...,N. 

The value of Mq n can be obtained from the requirement that the overlap coefficients 
provide a unitary transformation between the two bases. Using the orthogonality 
relation ( 1A.4I) of the dual —1 Hahn polynomials, we obtain 



(qi\m, N — m) 



Q m (x e ;a,/3;N), 



where Un-i is the weight function ( 1A.5I) of the dual —1 Hahn polynomials and N is an 
even integer. It is seen from the formula (IA.5I) of Appendix A that if fi x > —1/2 and 
fly > —1/2, the weight function u^-i is positive for all i G {0, . . . ,N}. The overlap 
coefficients obey the orthonormality relation 



N 



22{q e \m, N - m)(n,N - n\q e ) = 5 n 



1=0 



It now possible to obtain the overlap coefficients between the Cartesian and polar 
wavefunctions of the Dunkl oscillator. We first observe that the eigenstates of Q have 
the expansion 



l<?2n+p) 



l + (-l) P Cn+ P 

2 



\k,n + p; ++) + 



i-(-i) p C 



2 1 



\k,n + p; — ), 



with p = 0,1; the formula is also valid for n = p = 0. The inverse relations have the 
explicit form 



\k, n; ++) = 


[Cn-1" 

. Cn . 


|?2n-l) - 


"Cn + l" 
. 2zCn . 


l<?2n) 


\k,n; — ) = 


"Cn + l" 
L Cn . 


\<l2n-l) + 


"1-Cn" 

2z<^ n 


\0Hn) 



for These formulas can be used directly to obtain the overlap coefficients 

(k, n; + + \m, N — m), (k, n; \m, N — m), 

as linear combinations of dual —1 Hahn polynomials. 



4-2. Overlap coefficients for s x s y = — 1 

The overlap coefficients in the parity sector s x s y = —1 are obtained similarly to the case 
s x s y = 1. We again start by writing the energy eigenstates in the polar basis in terms of 



The Dunkl oscillator in the plane I 17 

the eigenstates of the operator J 2 . It follows from the relation (TIT?]) that the eigenstates 

of the operator J 2 with eigenvalues that we denote by |n, — h) j 2 and |n, ) j 2 are 

given by 

\ n ,-+)j 2 = a n = 2^(n + fi x )(n + fi y ), 

\n, ) j2 = -j=(\k,n;-+) + i|As, n; H — >), er^ = -2yJ(n + fj, x )(n + fj, y ), 

where n G {1/2, 3/2, . . .}. In this sector, the operator Q is equivalent to 

Q = iJ*Rx + (X ~ »y)Rx + (1/2)1 



The eigenstates \n, — ±)q and eigenvalues of Q are easily found to be 

1 

71 



n, -+)q = -i (£ n \n, -+)j 2 + \n, — )j 2 ) , g+ = 1/2 - 2n - /i x - /i y , 



and 

I n \ ^ = 

V2 



n, ) Q = —= (-Cn\n, -+)j 2 + \n, )j 2 ) , q n = 1/2 + 2n + fi x + ft y , 



where 

fi x - fi y + 2iy/(n + fi x ){n + ft y ) 
2n + fi x + fi y 

It is easily verified that £ n £* = 1. The eigenstates of Q can be grouped in a single 
expression. We write 

q e = {-l) t+1 (£ + ti x + fi v + l/2). 
The eigenvectors have the expressions 

kWs = ^((-i) p W |J + 1/2, -+}j 2 + |j + 1/2, —) j2 

We now compute the overlap coefficients between the eigenstates of Q and the 
eigenstates of "H expressed in the Cartesian basis for a given energy level £n- In the 
sector s x s y = — 1, the energy takes the values 

S N = N + fl x + fly + 1, 

where N is an odd integer. The eigenspace corresponding to Sn is spanned by the 
vectors 

|0,iV), |l,iV-l), ■■■\m,N -m), ■■■ , \N,0). 
We denote the overlap coefficients by 

{q e \m,N-m) = W* N . 
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The coefficients W l m N can be computed from the relation 

qi (qe\m, N — m) = (qi\ Q\m, N — m). (29) 
In terms of the parabosonic operators, the operator Q acting on the sector s x s y reads 
Q = {A X A\ - AlA y )R x + - iiy)Rx + (1/2)1 (30) 



Upon substituting (1301) in f l29l) and using the actions (1221) . one finds the recurrence 
relation 

^L,n = A m+ i W m+lfN + B m W l m N + A m W m _ 1>N , 

where 

A m = (-l) m ^[ml x [N-m + ll y , B m = (-l) m (m- - A*») + 1/2- 
After writing 

with Vo(qe) = 1 and introducing the monic polynomials 

V m {qe) = - A t-, 

■ ■ ■ Si m 

one finds that the polynomials V m {fLl) satisfy the three-term recurrence relation 

qeV m (qe) = Vm+i(qe) + B m P m (qe) + U n Pm-i(qt), (31) 

with 

U m = [m] ltt [N-m+ 

By comparing the recurrence relation (13T|) with that of the dual —1 Hahn polynomials 
(]A.3j) . one obtains 

V m (qt) = 2- m Q m (x e ,a,(3,N), 
with the parameter identification 

a = 2fi x , f3 = 2(i y , 

and the variable 

The requirement that the overlap coefficients provide a unitary transformation leads to 
the relation 

(qt\m, N-m) = J Wt Q m (xt, a, P, N), 

V U l ■ " " U m 

where N is an odd integer. The overlap coefficients satisfy the orthogonality relation 

JV 



y]{qi\m, N - m)(qi\n,N - n) = 5 nm . 



1=0 
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It is again possible to recover the overlap coefficients between the wavefunctions by 
expressing the eigenvectors of Q in terms of the eigenstates in the polar basis. One has 



l + (-l) p gj + l/2 
2 



|A;,j + l/2;-+) + 



(-l)%- + l/2-l 

2i 



\k,n;+-), 



with p G {0, 1}. The inverse relation reads 

|fc,j + l/2;-+>= [^f] \q 2j } + 
\kJ + l/2 ]+ -)= \tm±] \q 2j ) + 



2 *Sj + l/2 
2*^ + 1/2 



Wj+l), 

\Q2j+l)- 



Hence it is seen that the expansion coefficients between the Cartesian and polar bases are 
given in terms of linear combinations of dual —1 Hahn polynomials. These coefficients 
can also be expressed in integral form using the separated wavefunctions obtained in 
Section 2. 



5. The Schwinger-Dunkl algebra and the Clebsch-Gordan problem 

The Schwinger-Dunkl algebra and the dual —1 Hahn polynomials have both appeared 
in the examination of the Clebsch-Gordan problem for the Hopf algebra s/_i(2) [9j HO]. 
In this Section, we explain the relationship between the two contexts. This will clarify 
the introduction of the operator Q in the previous Section. 

5.1. sZ_i(2) Clebsch-Gordan coefficients and overlap coefficients 

The sZ_i(2) algebra is generated by the elements Aq, A± and R with the defining 
relations 

[A) j R] = 0, [A , At] = ±A±, {At, R} = 0, {At, A_} = 2A , 

and R 2 = I. It admits the Casimir operator 

Q = A + A^R - A R + (1/2)/?, 

which commutes will all the generators. This algebra has infinite-dimensional irreducible 
modules V^'^ spanned by the basis vectors i4 , n6E The action of the generators 
on the basis vectors is 

A v^ = (n + ii + 1/2)^, Rvj** = e(-l) n v^\ 

It is easily seen that Qvn'^ = —e/j,Vn'^. 

The s/_i(2) algebra is a Hopf algebra and has a non-trivial co-product. Upon taking 
the tensor product of two irreducible modules V^ 61 '^ 1 ' ® V^ 2 '^ spanned by the basis 
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vectors en <8> em 2 '^ 2 \ one obtains a third module V (in general not irreducible) by 
adjoining the action 

A (v ®w) = (A v) ®w + v<g) (A w), R(v ®w) = (Rv) ® (Rw), 

A±(v ®w) = (A±v) ® (Rw) + v® (A±w), 

where v G V^ 1 '^ and w £ V^ 2,fl2 \ On V, we have the Casimir element 

Q = (A^Af - A<pA!?>)BF) - (1/2)rVrW - e^R^ - e 2f i 2 R^ , 

where the superscripts indicate on which module the generators act; e.g. A± = A± (8)1. 

The eigenvalues of Q represent the irreducible modules y( e *>'*») appearing in the 
decomposition of V" = ® i V( eiW '. The Clebsch-Gordan coefficients of sZ_i(2) are the 
expansion coefficients between the direct product basis en <8>e&[ 2 and the eigenbasis 
/fc £l of the operator Q; this corresponds to the 'coupled' basis. Given the addition 
rule of A , one has 

ni+n2=A r 

where C^^In are the Clebsch-Gordan coefficients, which were shown to be given in 
terms of dual —1 Hahn polynomials in |U HO] . 

In our model, it is seen that the operators {H x , A x , A> x }, and {7i y , A y , A^ y } realize the 
two sZ_i(2) modules V 4 * and V 4 ", with e x = = 1. The Cartesian basis states |n 2 ,n a ) 
correspond to the direct product basis and the operator Q given in (1231) corresponds to 
the Casimir operator Q. This explains the origin of the operator Q in our approach to 
the overlap coefficients. 

5.2. Occurrence of the Schwinger- Dunkl algebra 

In our model, the Schwinger- Dunkl algebra occurs as the symmetry algebra. The algebra 
sd(2) also appears in the C.G. problem of s/_i(2) as a 'hidden' algebra. We illustrate 
how this comes about. 

In the C.G. problem, it follows from (132j) that the following operators act as multiple 
of the identity: 

A«+4 2) , Q {1 \ Q {2 \ R {1) R {2) - (33) 

In the direct product basis, in addition to the operators ( 133|) . the operators 

K = - 4 2) )/2, R = R {1) 

and R^ are also diagonal. In the 'coupled' basis, in addition to ( |33|) . we have the Casimir 
operator K\ — Q which is diagonal. Hence, the tensor product basis corresponds 
to having the operators (133]) plus K Q and R in diagonal form and the coupled basis 
corresponds to having the operators ( )33|) and K\ in diagonal form. A direct computation 
shows that the set {Kq, K\, R} generates the Schwinger-Dunkl algebra |9]. 

We have thus established the connection between our model and the Clebsch- 
Gordan problem of the algebra sl-\(2). 
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6. Conclusion 

We considered the Dunkl oscillator model and showed that it is a superintegrable system. 
We have exhibited the symmetry algebra that we called the Schwinger-Dunkl algebra 
and we have obtained the exact solutions of the Schrodinger equation in terms of Jacobi, 
Laguerre and generalized Hermite polynomials in Cartesian and polar coordinates. The 
expansion coefficients between the Cartesian and polar bases have been obtained exactly 
in terms of linear combinations of dual —1 Hahn polynomials and we established the 
connection between these overlap coefficients and the Clebsch-Gordan problem of the 
algebra s/_i(2). 

The representations of the symmetry algebra of a superintegrable system explain 
how the degenerate eigenstates of this system are transformed into each other. In the 
second series of the paper, we shall consider the representations of the Schwinger-Dunkl 
algebra. As will be seen, these representations exhibit remarkable occurrences of other 
— 1 polynomials. 

It would be of interest to consider in a future study the 3D Dunkl oscillator model, 
which will provide another example of a superintegrable system with reflections. It was 
shown in (TO] that the Bannai-Ito polynomials occur as Racah coefficients of the algebra 
sL_i(2). Given the connection between the 2D Dunkl oscillator and the Clecbsch-Gordan 
problem of sZ_j_(2), one can expect that the Bannai-Ito polynomials will occur in the 
description of the 3D Dunkl oscillator model. 
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Appendix A 

Appendix A.l. Formulas for Laguerre polynomials 
The Laguerre polynomials L^\x) are defined by [23J: 




(a + 1) 
n\ 



71 



a + l 



—n 
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where (a) n = (a) (a + 1) • • • (a + n — 1) is the Pochhammer symbol. They obey the 
orthogonality relation: 

e-VL^(x)L^(x) = r(n + ^ + 1) S nm , (A.l) 

n! 

for a > — 1. 



o 



Appendix A. 2. Formulas for Jacobi polynomials 
The Jacobi polynomials Pn^\x) are defined by [23J: 

-nn + a + (3 + l 1 — x 



a + 1 



p(^) (x)= (^±Ik 2Fl 
n! 

They obey the orthogonality relation: 

[ n - rfn 4- t^P^VtIP^M - 2 ° +/?+1 r(n+a+l)r(n+/3+l) r /a o\ 

y ^ A X / V L ^ X ) r m \- h ) r n \ x ) — 2n+a+/3+l r(n+a+/3+l)n! l^-^J 

provided that a > — 1 and /3 > — 1. 

Appendix A.3. Formulas for dual — 1 Hahn polynomials 

The monic dual —1 Hahn polynomials Q n (x; a, (3; N) have the recurrence relation [41j: 

xQ n {x) = Q n+1 {x) + b n Q n (x) +u n Q n _i{x). 
The recurrence coefficients are given by: 

At i r 7i r X11 , /(-l)"' +1 (2e + 2C)-l, iVeven, 

u n = 4 \n \e \N — n + 1 U b n = < , (A.3) 

L J?L K \(-l)"(2C-20-l, AT odd, 1 ' 

where 

_J^f^i, Neven, _ N even, 

~{f, A^odd, ' ^ \§, A^ odd. 

They obey the orthogonality relation: 

N 

S ^2uj t Q n {xt)Q m {x ( >) = v n 5 nm (A.4) 
t 

The weight is given by 

{ (-\)i{-m) ]+q (l-a/2)j (1-^/2-/3/2), (l-/3/2) m AT p V pn 

j! (1-/3/2), (m+l- a /2-/3/2) J+9 (l-a/2-/5/2) m ' iveven, 

(-l) J( - m)j (l/2+a/2) J+q (l/2+a/2+/3/2) J (l/2+/3/2) m+1/2 
j! (1/2+^/2) j+ ,(m+3/2+a/2+ / 9/2) J - (l+a/2+/3/2) m+1/2 > UUU ' 

where g G {0, 1} and with m = N/2, v n — U\ ■ ■ ■ u n . The grid points have the expression 

(-lY(2£ + l-a- /3), Neven, 
-l) e (2e + l + a + /3), A^ odd. 
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We here indicate how it can be simply seen that the Dunkl derivative (j2J) is anti- 
Hermitian with respect to the scalar product (jSJ). This is recorded for completeness 
and convenience. We need to check that 



Tl)l{x)[V%i) X {x)]\x\ 2 ^dx 



[V^ 2 (x)]*^M 2flx dx = -(V^M, 



for all functions ipi(x), ip 2 {x) belonging to the L 2 space associated to the scalar produt 
([H]). We split the computation in the four possible parity cases for ipi(x), ip 2 {x). This is 
sufficient since any function can be decomposed into its even and odd parts and since 
the scalar product §5§ is linear in its arguments. 

In the even-even case, one has ipi(x) = tpi(-x), ip 2 (x) = ip 2 (—x). It follows that 



r 2 (x)[D^ l (x)}\x\^dx= / ^(x)[d x Mx)\\x\ 2>1 dx = 0, 



since the integrand in odd. Similarly, we have {TJ^ip^ipi) = 0. 

In the odd-odd case ip\{x) = —ipi(—x), ip 2 {x) — —^(—x) and one obtains 



ip^p^x^Wx^dx 



d x ipi{x) + —ij)i(x) 
x 



:\ 2fl dx = 0, 



since the integrand is odd. Similarly, we have (TJ^^lipi) = 0- 

In the even-odd case, ipi{ x ) = ~i , i{~ x )i ^2 0*0 — 4>2{—x) and it follows that 



r 2 (x)[v^Mx)]\x\ 2fi dx 



2a 

dx4>i(x) H ipi(x) 

x 



\x\ 2fl dx 



r 2 (x) 



2u 

dxipi(x) H ipi(x) 

x 



2^i{x)^)* 2 {x)x 2 ^ 



r 2 (x) 

x 2[1 dx 

d x [i)\{x)x 2 ^ x {x)dx + 4a I ^{x^^x^^dx 



[d x r 2 (x)]Mx)\x\ 2ll dx 



[V^ 2 (x)]*Mx)\x\ 2f "dx, 



where we have used the vanishing conditions on ipi{x), ip 2 {x) at infinity. It thus seen 

that (V> 2 |Z>£Vi> = -(VM^i). 

In the even-odd case, one has ip\{x) = ipi(-x), ip 2 (x) = —^(—x) an d one obtains 



il>Z(x)[T>£Mx)]\x\ 2 »dx 



^ 2 (x)[d x ipi(x)]\x\ 2fM dx = 2 1 ^2(x)[d x tj}i(x)} x 2f *dx 



2ipi(x)if)2(x)x 



2fl 



8 ^ W + % ;W 

X 



ipi(x) x 2>1 dx 



[D^ 2 (x)]*Mx)\x\ 2tl dx, 



where we have used the vanishing conditions on ipi(x), ip 2 (x) at infinity. Hence we have 
(^2|£>^i> = -(V%ip 2 \ipi) and the result 



is established in all cases. 
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